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ABSTRACT 

The distribution of satellite galaxies relative to isolated host galaxies in the Sloan Digital Sky Survey 
(SDSS) is investigated. Host-satellite systems are selected using three different methods, yielding samples 
of ~ 3300, ~ 1600, and ~ 950 satellites. In the plane of the sky, the distributions of all three samples 
show highly significant deviations from circular symmetry (> 99.99%, > 99.99%, and 99.79% confidence 
levels, respectively), and the degree of anisotropy is a strong function of the projected radius, rp, at 
which the satellites are found. For rp <^ 100 kpc, the SDSS satellites are aligned preferentially with 
the major axes of the hosts. This is in stark contrast to the Holmberg effect, in which satellites are 
aligned with the minor axes of host galaxies. The degree of anisotropy in the distribution of the SDSS 
satellites decreases with rp and is consistent with an isotropic distribution at of order the l-a level for 
250 kpc ^ r-p < 500 kpc. 

Subject headings: galaxies: fundamental parameters — galaxies: halos — galaxies: structure — dark 
matter 



1. INTRODUCTION 

The distribution of satellite galaxies relative to their 
host galaxies may hold important clues to the formation 
history of large galaxies. It should reflect the time scales 
for dwarf galaxy evolution and orbital decay, as well the 
nature of the host's dark matter halo and the manner in 
which the host accretes mass. On scales of order the virial 
radius, r2oo, the halos of large galaxies in cold dark mat- 
ter (CDM) universes are flattened with a mean projected 
ellipticity of order 0.3 (e.g., Dubinski & Carlberg 1991; 
Warren et al. 1992; Jing & Suto 2002). On smaU scales, 
r << r2oo, gas cooling makes the halos somewhat rounder, 
but on scales r ~ r2oo baryons do not affect the shapes of 
the halos very much (e.g., Kazantzidis et al. 2004). Evi- 
dence for non-spherical halos includes dynamics of polar 
ring galaxies, geometry of X-ray isophotes, flaring of HI 
gas in spirals, evolution of gaseous warps, kinematics of 
Population II stars in our own Galaxy, and gravitational 
lensing (se, e.g., Sackett 1999; Keeton et al. 1998; Mailer 
et al. 2000; Hoekstra et al. 2004). However, Ibata et al. 
(2001) concluded a nearly spherical halo for the Milky Way 
on the basis of tidal stellar streams. 

If the halos of large galaxies are flattened, the distribu- 
tion of their satellites may reflect the flattening of the halo 
and, hence, show some preference for alignment with the 
major axis of the halo mass distribution, as projected on 
the sky. Assuming the light and mass of a large galaxy are 
aligned reasonably well, one then expects a preference for 
satellites to be aligned with the major axes of the images 
of their hosts. However, this has not generally been ob- 
served. Holmberg (1969) investigated the distribution of 
satellites relative to disk galaxy hosts and, for projected 
radii rp <^ 50 kpc, found a tendency for the satellites to 
be located preferentially close to the minor axes of the 
hosts. In addition, the 11 brightest satellites of the Milky 
Way form a flattened structure that is oriented roughly 



perpendicular to the disk (e.g., Lynden-Bell 1982; Ma- 
jewski 1994). Conversely, Valtonen et al. (1978) found a 
tendency for compact satellites to be aligned with the ma- 
jor axes of disk hosts. Hawley & Peebles (1975), Sharp et 
al. (1979), and MacGillivray et al. (1982), suggested that 
any anisotropy in the distribution of satellite galaxies was 
at most rather small, and perhaps non-existent. Zaritsky 
et al. (1997), hereafter ZSFW, studied the distribution of 
satellites around spiral hosts and found evidence for an 
alignment of the satellites along the host minor axes for 
300 kpc ^rp 500 kpc. Sales & Lambas (2004), hereafter 
SL, used a set of 1498 host galaxies with 3079 satellites 
from the Two Degree Field Galaxy Redshift Survey (2dF- 
GRS; CoUess et al. 2001, 2003) to investigate the distribu- 
tion of satellites and found a preference for the satellites 
to be aligned with the minor axes of the hosts, but only 
for a very specific subsample of host-satellite pairs. 

The theory regarding the Holmberg effect is equally as 
murky as the observations. ZSFW were unable to recover 
their observed anisotropy from CDM simulations. Abadi 
et al. (2003) suggested that the Holmberg effect could be 
caused by the cumulative effects of accretion of satellites 
by the host. Peharrubia et al. (2002) suggested that the 
Holmberg effect orginates with the fact that polar orbits 
inside a massive, flattened halo do not decay as quickly as 
planar orbits. However, Knebe et al. (2004) found that the 
orbits of satellites of primary galaxies in cluster environ- 
ments were located preferentially within a cone of opening 
angle 40° (i.e., planar alignment, not polar). Since the 
structure of CDM halos is largely independent of the halo 
mass scale (e.g., Moore et al. 1999), the implication of this 
result is a preference for the satellites of isolated galaxies 
to be aligned with the major axes of the host halos. In- 
deed, the most recent numerical work seems to show that 
this the case (e.g., Agustsson & Brainerd 2005, hereafter 
AB; Libeskind et al. 2005; Zentner et al. 2005). 
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The organization of this Letter is as follows. The distri- 
bution of satellite galaxies relative to host galaxies in the 
third data release (DR3) of the SDSS (Abazajian et al. 
2005; Pukugita et al. 1996; Hogg et al. 2001; Smith et al. 
2002; Strauss et al. 2002; York et al. 2002) is shown in §2 
and is compared to an isotropic distribution. The results 
obtained in §2 arc compared to other recent observational 
and numerical results in §3. Throughout, the null hypothe- 
sis is that satellites are spherically distributed about their 
hosts, giving rise to a circularly-symmetric distribution 
on the sky. A rotation of a spherically symmetric dis- 
tribution through any combination of Eulcr angles always 
gives rise to a 2-d distribution that is circularly-symmetric 
and, hence, any deviation in the satellite distribution from 
pure circular symmetry cannot be caused simply by pro- 
jection and/or rotation effects. Cosmological parameters 
Oo = 0.3, Ao = 0.7, and Hq = 70 km sec~^ Mpc ~^ are 
adopted below. 

2. HOST-SATELLITE SELECTION AND THE DISTRIBUTION 
OF SATELLITES 

The r' band photometric data from the SDSS DR3 

were used to define the apparent magnitudes of the galax- 
ies, their position angles, and the shapes of the equiv- 
alent image ellipses. Only galaxies for which the SDSS 
redshift confidence parameter, zconf, is greater than 0.9 
were used. Three methods were adopted in order to se- 
lect isolated host galaxies and their satellites. "Sample 
1" was selected using the criteria of SL: [1] Hosts are at 
least 2.5 times brighter than any other galaxy within a 
projected radius of rp < 700 kpc and a relative radial 
velocity difference of \dv\ < 1000 km sec^^. [2] Satel- 
lites are at least 6.25 times fainter than their host, are 
found within Vp < 500 kpc, and the host-satellite veloc- 
ity difference is \dv\ < 500 km sec~^. "Sample 2" was 
selected using the criteria of McKay et al. (2002) in their 
study of the dynamics of satellite galaxies in the SDSS, 
with the exception that here a slightly smaller maximum 
Tp for the satellites was used in order to be consistent 
with the other two samples: [1] Hosts are at least 2 times 
brighter than any other galaxy within Vp < 2.86 Mpc and 
\dv\ < 1000 km s~^. [2] Satellites are at least 4 times 
fainter than their host, are found within Vp < 500 kpc 
(cf. rp < 714 kpc in McKay et al. 2002), and have 
|fiw| < 1000 km s~^. "Sample 3" was selected using the 
criteria of Zaritsky ct al. (1993) and ZSFW: [1] Hosts are 
at least 8 times brighter than any other galaxy within 
rp < 500 kpc and \dv\ < 1000 km sec~^. In addition, 
hosts are at least 2 times brighter than any other galaxy 
within rp < 1 Mpc and \dv\ < 1000 km sec~^ [2] Satel- 
lites arc at least 8 times fainter than their host, are found 
within rp < 500 kpc, and have \dv\ < 500 km sec~^. 

Additionally, in order to eliminate hosts which have a 
very large number of satellites (and may be associated 
with cluster systems), the luminosity of a host must be 
greater than the sum of the luminosities of its satellites. 
Finally, the samples were restricted to hosts with ellipticity 
is e > 0.2 so that the position angles of the major axes of 
the hosts arc well determined. The selection criteria yield 
1973 hosts and 3292 satellites for sample 1, 913 hosts and 
1575 satellites for sample 2, and 575 hosts and 935 satel- 
lites for sample 3. The median redshifts of the hosts are 
Zmed = 0.05, 0.08, and 0.03, respectively. The number of 



satellites per host is shown in Fig. 1, from which it is clear 
that the samples are dominated by systems containing one 
or two satellites. Also, although the morphologies of all of 
the hosts are not known, a visual inspection of images of 
the 50 brightest hosts (all of which appear in all 3 samples) 
shows that 44 are spirals (mostly Sc or Sd), 4 are SO, and 
2 are Sd/Irr. 
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Figure 1. Normalized probability distribution of the number of satellites 
per host for each of the three samples of hosts and satellites (see text for 
sample selection criteria). 

The distribution of the satellites relative to their hosts 
was computed by measuring the angle, (/), between the ma- 
jor axis of the host and the direction vector on the sky that 
connects the centroids of the host and its satellite. Because 
we are simply interested in investigating any preferential 
alignment of the satellites with the principle axes of the 
hosts, (f) is restricted to the range [0°,90°] (i.e., (f> ^ 0° 
indicates alignment with the host major axis and </> = 90° 
indicates alignment with the host minor axis). Results of 
the analysis for all three samples are shown in Fig. 2. The 
left panels show the differential probability distribution, 
P((/)), and the right panels show the cumulative probabil- 
ity distribution, P{(j) < 0MAx)- The probabilities expected 
for an isotropic (i.e., circularly-symmetric) distribution of 
satellites are shown by the dotted lines in all eight panels of 
Fig. 2. A test was applied to P{4>), and an isotropic dis- 
tribution of satellites was rejected for all three samples (see 
Fig. 2). The mean value of 0, (0), was also computed for 
each sample of satellites and excellent agreement amongst 
the samples is found. The mean value of (j) is considerably 
less than the expected value for an isotropic distribution 
(i.e., {(j)) = 45°), and differs formally from the isotropic 
distribution by 6.4-cr, 4.6-cr, and 3.3-cr for samples 1, 2, 
and 3 respectively. Finally, a Kolmogorov-Smirnov (KS) 
test was applied to P{(j) < ^max) for each of the sam- 
ples and, again, in all cases an isotropic distribution was 
rejected with high significance (see Fig. 2). 

Fig. 3 shows the dependence of {<p) on the projected 
separation between the satellites and their hosts. For a 
given sample, the widths of the bins were chosen such that 
an equal number of satellites fell into each bin. Again, 
the dotted line indicates the expectation of an isotropic 
distribution. From Fig. 3 it is clear that the anisotropy 
in the distribution of the satellites is strongest on rel- 
atively small scales {vp ^ 100 kpc), and on scales of 
250 kpc ^ Tp ^ 500 kpc all three samples are consistent 
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with a uniform distribution at of order the l-a level. It 
is very unlikely that on scales Vp <J 100 kpc a large num- 
ber of giant H-II regions (i.e., analagous to 30 Doradus) 
have been misidentified as "satellites" for a number of rea- 
sons. First, only objects for which the SDSS spectra are 
clearly classified as "galaxies" have been used. Also, the 
ratio of the luminosity of the hosts to the luminosity to 
these objects is much smaller than would be expected for 
giant H-II regions (median luminosity ratios of 13, 10, 
and 24 for samples 1,2, and 3 respectively). Finally, the 
isophotal areas of these objects compared to their hosts is 
much larger than would be expected for giant H-II regions 
(median isophotal area ratios of 0.09, 0.12, and 0.04 for 
samples 1, 2, and 3 respectively). 
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Figure 2. Left panels; Differential probability distribution for the loca- 
tion of satellites relative to the host major axis. Error bars were eom- 
puted from 1000 bootstrap resamplings of the data. Dotted line shows 
P(0) for an isotropic distribution. Formal rejection confidence levels 
from the test are shown in the panels, along with {(j)) for each of 
the samples. Right panels: Cumulative probability distribution for the 
location of satellites relative to the host major axis. Formal rejection 
confidence levels from the KS test are shown in the panels. 

3. DISCUSSION AND CONCLUSIONS 

All three samples of satellites show a strong preference 
for being aligned with the host major axis on small scales 
{rp <^ 100 kpc), while on larger scales (250 kpc <, rp <^ 
500 kpc) the satellite distribution is fairly consistent with 
an isotropic distribution. These results contrast sharply 
with most previous studies; i.e., when an anisotropy in 
the distribution of satellites has been seen in the past, the 
satellites have generally been aligned with the host mi- 
nor axis, not the host major axis. Aside from this Letter, 
the most recent investigation of the distribution of satel- 



lite galaxies relative to host galaxies is that of SL, who 
concluded that the satellites of 2dFGRS host galaxies arc 
preferentially aligned with the minor axes of the hosts, 
but only for a very select subset of the systems. That is, 
SL selected hosts and satellites in exactly the same way 
in which sample 1 of this Letter was selected, and host- 
satellite pairs for which the relative velocity difference is 
\dv\ < 500 km sec~^ were specifically included in SL's 
work. When SL examined their entire sample of satellites, 
no anisotropy in the distribution of satellites was found. 
When they restricted the sample to hosts and satellites 
with \dv\ < 160 km sec~^, a strong detection of the Holm- 
berg effect was seen. In particular, SL report that the 
ratio of the number of satellites found within 30° of the 
host major axis to the number of satellites found within 
30° of the host minor axis is / = 0.80 ± 0.04 for hosts and 
satellites that have low relative velocity. In contrast, the 
data presented here yield / = 1.30 ± 0.08 for all of the 
satellites in sample 1 and / = 1.28 ± 0.10 for hosts and 
satellites in sample 1 for which < 160 km sec~^. 
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Figure 3. Mean orientation of satellites with respect to the host major 
axis as a function of the projected radius. 

The SL result is troubling because the typical error on 
a radial velocity measurement in the 2dFGRS is of order 
85 km sec~^, so a measured value of \dv\ = 160 km sec~^ 
is comparable to the error in the value of \dv\. (In the case 
of the SDSS, the typical velocity errors are of order 20 to 
30 km sec"""^.) In addition, simulations of galaxy redshift 
surveys by van den Bosch et al. (2004) show that the frac- 
tion of interloper galaxies is higher for host-satellite pairs 
with small values of \dv\ than it is for host-satellite pairs 
with large values of Id?;!. It is possible, therefore, that the 
host-satellite pairs in SL with \dv\ < 160 km sec~^ are 
heavily contaminated by interlopers and some unknown 
selection bias is responsible for the signal which they de- 
tect. Further, the results in this Letter, combined with the 
results of van den Bosch et al. (2005), may provide some 
insight into the apparent lack of anisotropy seen by SL in 



4 



Anisotropic Distribution of SDSS Satellite Galaxies 



their full sample. Fig. 3 of this Letter indicates that the 
anisotropy occurs primarily on small physical scales. How- 
ever, when van den Bosch et al. (2005) combined mock 
galaxy redshift surveys with the 2dFGRS, they found a 
marked absence of satellites at small projected radii, which 
they attributed in part to the merging and overlap of im- 
ages in the APM catalog. That is, SL might have had too 
few host-satellite pairs with ^ 100 kpc to detect the 
anisotropy that is seen here in the SDSS DR3 samples. 

Following on from the numerical work of Knebe et al. 
(2004), who first predicted a preferential alignment of 
satellites with the major axes of primary galaxies in clus- 
ter environments, the most recent theoretical work on the 
distribution of the satellites of isolated host galaxies shows 
that in a CDM universe the satellites ought to be found 
preferentially close to the principle plane of the host halo 
(i.e., the plane defined by the two largest moments of in- 
ertia of the host halo). High-resolution simulations of the 
formation of "Milky Way" systems show that the satel- 
lite galaxies form a rather thin, planar structure that lies 
within the principle plane of the halo (e.g., Libeskind et al. 
2005; Zentner et al. 2005). In addition, AB showed that, in 
projection on the sky, when hosts and satellites are selected 
from a ACDM simulation that incorporates semi-analytic 
galaxy formation, the satellites show a strong preference 
for clustering near to the major axes of the projected host 
halos. 

AB placed circiilar disks in the principle planes of the 
host halos (i.e., to represent the disks of the host galax- 
ies since the simulation did not contain images of visi- 
ble galaxies) and they then computed the distribution of 
satellites with respect to the major axes of the projected 
circular disks using 100 random rotations of the simula- 
tion. AB found that, in the plane of the sky, satellites 
with rp < 500 kpc showed a clear preference for cluster- 
ing nearby to the major axes of the projected circular disks 
and, in fact, the degree of anisotropy shown by AB's satel- 
lites was somewhat greater than the degree of anisotropy 
shown by the SDSS satellites. The anisotropy shown by 
AB's satellites is caused by the fact that in a CDM uni- 
verse satellites are accreted preferentially into the principle 
plane of the host halo along filaments (see also Libeskind 
et al. 2005 and Zentner et al. 2005). 

AB were, however, unable to reproduce the trends 
shown in Fig. 3; i.e., that the anisotropy in the distribution 
of the SDSS satellites arises primarily on scales 100 kpc. 
Contrary to the results of this Letter, AB found that the 
anisotropy in the distribution of the simulated satellites 
persisted to large scales (r^ ~ 500 kpc). Since the hosts 
and satellites in the simulation were selected in exactly the 
same way as was done here for the SDSS, the fraction of 



galaxies in the SDSS that may have been falsely identified 
as satellites (i.e., "interloper galaxies" ) should be similar in 
both cases. Hence, the well-known increase in the number 
of interlopers with increasing Vp (e.g., Prada et al. 2003; 
Brainerd 2004) is unlikely to be the cause of the SDSS 
satellites showing a nearly isotropic distribution at large 
Tp. It is not clear what the underlying cause of the discrep- 
ancy between the degree of large-scale anisotropy shown 
by AB's satellites and the SDSS satellites happens to be, 
but a likely possibility is that the disks of host galaxies 
may be aligned only with the very innermost parts of their 
halos (whereas AB used all mass within the virial radius 
to define the plane in which the host's disk would reside). 
Recent work by Bailin et al. (2005) on the formation of 
disk galaxies in a CDM universe shows that although the 
orientation of the disk is aligned well with the halo on 
small scales, the orientation of the disk is essentially un- 
correlated with the outer regions of the halo. Since the 
satellites presumably trace the gravitational potential and 
accretion history of the halo (i.e., the mass), rather than 
the disk specifically, Bailin et al.'s result would lead di- 
rectly to the expectation that on very large physical scales 
the satellites of disk galaxies should be distributed roughly 
isotropically relative to the direction vector defined by the 
major axis of the disk. 

In conclusion, the satellites of host galaxies in the SDSS 
DR3 show a marked preference for alignment with the ma- 
jor axes of the hosts. Although the sense of this anisotropic 
distribution agrees well with the expectations of a CDM 
universe, a proper comparison of the details of the signal 
(i.e., a strongly anisotropic distribution on small scales and 
a nearly isotropic distribution on large scales) awaits high 
resolution simulations that are capable of fully addressing 
the orientations of visible galaxies within their dark matter 
halos. 
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